Spin-Hall effect in a disordered 2D electron-system 
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We calculate the spin-Hall conductivity for a two-dimensional electron gas within the self- 
consistent Born approximation, varying the strength and type of disorder. In the weak disorder 
limit we find both analytically and numerically a vanishing spin-Hall conductivity even when we 
allow a momentum dependent scattering. Separating the reactive from the disspative current re- 
sponse, we find the universal value afjj — e/Svr for the reactive response, which cancels however 
with the dissipative part crf^ = — e/Syr. 



Spin-orbit coupling in two dimensional electron sys- 
tems allows a number of unconventional transport phe- 
nomena, since charge current and the spin degrees 
of freedom are coupledi. In particular the spin- 
Hall effect in two-dimensional electron systems, i.e. 
a spin-current which flows in the plane but perpen- 
dicular to the electrical current, and which is polar- 
ized perpendicular to the plane, has been discussed 
intensively^i3,4,5,6,7,8,9,io,ii,i2,i3,i4,i5,i6,i7 o^gj. tlje 1^^^^ 

year. The spin-Hall conductivity connects the spin- 
current with an electric field, jy = asnEx, where jy 
denotes a current in y-direction with spin-polarization 
in z-direction. In a clean two dimensional electron gas, 
the spin-Hall conductivity was predicted to have a uni- 
versal value — g/Stt, independent of the strength 
of the spin-orbit scattering^. Several publications have 
addressed the issue of whether this result is modified in 
the presence of impurity scattering. Murakam:^^ ana- 
lyzed the Luttinger-Hamiltoniani^, which applies to two- 
dimensional hole gases, and concluded that the spin-Hall 
conductivity in the limit of weak impurity scattering re- 
produces the intrinsic value (at least when restricting to 
s-wave impurity scattering). For the Rashba modeliS, 
which applies to two-dimensional electron gases, con- 
flicting results exist in the literature: By applying the 
standard Green's function techniques Inoue et al^ and 
Mishchenko et al>i^ concluded that s-wave impurities 
suppress the spin-Hall effect in bulk samples even when 
the disorder broadening of the energy levels is small com- 
pared to the spin-orbit splittingii. On the other hand, 
Dimitrova— and Chalaev et alii^, starting from the same 
model Hamiltonian and applying similar methods, found 
a non-zero spin-Hall conductivity. Even direct numeri- 
cal evaluations of the effect do not fully agree with each 
other: Xiong and Xie^ found within a scattering matrix 
approach the universal value of the spin-Hall conductance 
GsH = e/87r over a large parameter range. Nomura et 
aliiS on the other hand found a spin-Hall conductivity of 
the order of but not identical to the universal value. 

In this paper we calculate the spin-Hall conductivity 
for a bulk sample within the self-consistent Born approx- 



imation. We confirm Refs&iisi^, i.e., we find that even a 
weak disorder suppresses the spin-Hall conductivity. For 
s-wave scatterers we calculate the impurity self-energy 
and the dressed current vertex numerically. This allows 
us to obtain results beyond the limit epT — > 00, which is 
accessible analytically. We find that a non-zero spin-Hall 
conductivity is, in principle, possible although it remains 
much smaller than c/Stt. 

Our calculations are based on our previous workSSiSi, 
where a number of technical details can be found. In 
the following we sketch the derivation of the spin-Hall 
conductivity. The starting point is the Hamiltonian 



p 

H = h aa • p X 
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where the parameter a describes the strength of the spin- 
orbit coupling, (7 is a vector of Pauli matrices, and is a 
unit vector perpendicular to the two dimensional system. 
The spin-Hall conductivity is obtained by the standard 
linear response theory as 
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with G<(e) = /(e)(G^-G^), /(e) being the Fermi func- 
tion. In Eq. (0), the spin- and charge-current operators 
are given by j^. = (1/4){(72V -I- vcr^} and jc = v, respec- 
tively. The velocity operator, v, is obtained from the 
Hamiltonian ^ and reads w^'^ = p^'^ /m =F aay^x- We 
choose the electron charge as — e (e > 0). The trace in 
Eq. ^ is over the eigenstates of the Hamiltonian, and 
the bar indicates that the expression must be averaged 
over the disorder configurations. 

When performing the disorder average, we rely on the 
self-consistent Born approximation. To begin with, we 
consider point-like, i.e., pure s-wave scatterers. The re- 
tarded/advanced impurity self-energy is then given by 
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G«'^(p). 
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Due to the spin-orbit coupling, the Green's functions 
have a non-trivial structure in the spin-space, although 
the self-energy remains diagonal. Explicitly, one finds 
that Sss' = So'^ss', Gss' = GqSss' +Gia^„, -f G2CT^,, with 



Go(p) = i(G+ + G_) 

Gi(p) = i^(G+-G_) 
2 p 

G2(P) = -i^(G+-G_) 
2 P 

/ p2 

G± = e -I- /X - — =F - So 
V 2m 



(4) 
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By taking the zero frequency limit of Eq. (0), the spin- 
Hall conductivity reads 

'^^^ = E [2jTG«(p) J,-G^(p)] , (8) 
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FIG. 1: The dressed vertex Ff in units of a as function 
of disorder strength l/(eFT). Ff enters the dressed charge 
current, = Px/m + T^cjy and thus the spin-Hall conduc- 
tivity, Eq. 1121 . In comparison with the bare chare current, 
= Px/m ~ auy, the spin-dependence is strongly reduced. 




since terms of the type G^G^ and G^G^ contribute 
only in the order (l/ei?r)(a/vF)^ and can be safely ne- 
glected in the limit app ^ ef and/or for weak disorder 
epT ^ 1. The charge current has to be calculated in- 
cluding the vertex corrections, = p^/m-t-T^, compare 
Eq. (33) of Ref.^^. In the case of s-wave impurity scatter- 
ing, the momentum dependent part of the current vertex 
is not renormalized, while the momentum independent, 
but spin-dependent part, F^, is obtained by solving the 
set of equations 



ss' ! ss' 
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with the efFective bare vertex given by 
1 



j;^, = -aa;^, + — — V G^ '-^^^ 
2ttNqt ^ 
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By expanding F^^, — T^cr^^, in Pauli matrices we 
obtain the spin-Hall conductivity as 



^F^ImV^Go«(p)Gf(p) 



p 



(11) 



Performing the momentum integration Eq. 111() under 
the restriction that oipp ^ ef and tpT 3> 1 leads to 

<^sH = --i2 7riVoT (12) 

TT 1 + da^pJnT^ 



i.e., the universal value cr^// — e/din is recovered in the 
weak disorder limit. On the other hand, by inserting 
the dressed vertex FS w as calculated in RefsiS242i one 



finds that asH 
found ii&iiJ4. 



0. Our result then agrees with that 



As explained in Rcf.^i, the vanishing of the dressed 
vertex Ff is due to the fact that the integral on the right- 
hand side of Eq. (fTUIl gives ~ aa^ , making the effective 
bare vertex 7 itself to vanish. A more careful numerical 
evaluation of the integral for arbitrary disorder strength 
actually shows that the compensation of the two terms 
in Eq. I|ll)fl is exact only in the weak disorder limit, i.e., 
eFT ^ 1. In Fig. ^ we show the dressed vertex as a 
function of disorder. F| goes to zero as l/epT 0, 
nothing special is observed as app ~ I/t", and even in 
the strong disorder limit Fj remains much smaller than 
its bare value {—a). We conclude that although, in prin- 
ciple, a non-zero spin-Hall conductivity may be obtained, 
one expects a much smaller value than the universal one. 

Next we address the question whether the spin-Hall ef- 
fect is sensitive to the type of disorder potential. Inoue 
et alJi argued that asH niay be non-zero for long-range 
defect potentials, although an explicit result has not been 
given. In our calculation we follow again closelySi where 
all the details can be found. Here we assume weak dis- 
order, so that the inequalities ep ^ app ^ 1/t hold. 
In the following, we work in the eigenstate basis of the 
Hamiltonian |^ 



Apparently goes to zero when the spin-splitting, 
app, is small compared to the disorder broadening of 
the levels, l/r. If we neglect vertex corrections, i.e., if 
we insert in Eq. Hll|l the bare vertex Ff = —a, we find 



|P±> 



V2 



{±iexp(-i(^)|pT) + |pi)} (14) 



e Aa^p'^pT^ 
8nl + 4a2p2 T-2 



bare vertex 



(13) 



where tan((/j) = Py/px and the corresponding eigenvalues 
are E± = p'^ /2mzLap. In this basis, the matrix elements 



3 



of the current operators read 
(p±|j^|p±) - 
(p±|j-|pT) = 

(P±|jc|P±) 

(ptIj:Ip±) 



2 TO 

( — ± a] cos((i?) 
Vm / 

^ia sin((/?). 



(15) 
(16) 

(17) 
(18) 



To use Eq. ^ we need the dressed charge operator, J^. 
Since, as seen from the above equations, the spin-current 
operator is off-diagonal in the eigenstate basis we get the 
spin-Hall conductivity in the form 

'^sH = -^J2^^ [(P + bIIP-)(P - I Jc |P+)G^(P)G^(P) 

(19) 

To calculate the dressed current operator we make use of 
the assumption that the spin-orbit splitting is large com- 
pared to the impurity broadening of the levels, app ^ 
l/r. The off-diagonal matrix elements of the current op- 
erator are then obtained in terms of the diagonal ones, 

(pt|Jc1p±) = (ptIj.1p±) 



■E 



[{pT\V\p'm) (p'TO|y|p±) 



G^,ip')Gt{p'){p'm\J^\p'r 



(20) 



The diagonal matrix elements on the other hand were 
already considered in RefiSi, and are obtained from the 
equation 



(p±|J,-|p±) = (p±|j:|p±) 



(21) 



p',m 



\{p±\V\p'm)\'G^ip')Gi{p'){p'm\j:\p'm). 



We consider impurity scattering which conserves spin, 
but allow the scattering amplitude to be momentum- 
transfer dependent. Such a dependence appears as a 
product of two contributions. The first is due to the 
type of disorder potential one considers, V^,p', while the 
second is induced by the transformation to the eigenstate 
basis. The latter gives rise to the following matrix ele- 
ments 



(p - |y|p'±) (p' ± |y|p+) = - v')|Vp,pf (22) 



and 



\{p±\V\p'±)\' = -|yp^p,p(l + cos(^-^')) (23) 
|(p±|y|p'T)P - i|yp^p.p(l-cos(^-^')) (24) 

We assume that the scattering amplitude Vp_p' depends 
weakly on the momentum transfer, so that the scatter- 
ing probability depends only on the angle between the 



incoming and scattered particle. Under this condition 
we can expand the scattering probability as 

I Vp p' p = Vo + 2^1 cos(</7 - if') -f IV2 cos(2v9 - 2(^') + • ■ ■ ■ 

(25) 

To obtain the spin-Hall conductivity the two momentum 
integrations over p and p' have to be performed. We 
split the momentum integration in an integral over the 
energy ^ — /2m — fi and the angular variable (p, 



and find 



d£_G^G^ = 







(26) 



(27) 



2app + i/r app 
No J d^GfC^ = 2TrN±T±, (28) 

where N± and t± arc the density of states and the life- 
time in the two subbands. Notice that the first of the 
two integrals is correct only to the lowest order in the 
(small) parameter app/tp, whereas the second integra- 
tion is valid beyond that limit. Finally the spin-Hall 
conductivity is determined as 



e 

8^ 
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1 N+T+Jl - N^T^J-^ V2 - Vo 



4a 



Not 



Vo 



,(29) 



where the second term is due to vertex corrections. When 
calculating the product N±t±J^ to zero order in app /ep 
the vertex corrections disappear. Expanding the density 
of states, the scattering time and to dressed current op- 
erator to first order yields 



7V± 



No [IT 
r I 1 ± 



Vo 



Vo~Vi m 



app 
2ep 
V\ app 
Vo'2^ 



Pf Vo 
TO. 



Vo 



Vo~V2 



(30) 
(31) 
(32) 



Notice that the dressed current operator, to the leading 
order in a, is of the familiar form J — jt^c/t, where 
Ttr is the transport scattering time. By combining all 
the terms, one then finds that the spin-Hall conductivity 
(|29|l vanishes as in the case of pure s-wave scattering, 
<JsH = 0. 

As a last useful observation, we separate the reactive 
and dissipative contributions to the current response, 



sH 



where 



sH 



= lim — / — Tr 
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jfG<(e)j-ReG^(e-c.) 



+jyRcGR{e)]-G<{e-uj) 
e 



-Tr 



j|'ImG«j;fImG« 



(33) 
(34) 
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Since the zero frequency spin-Hall conductivity is real, 
the terms with imaginary (real) current matrix elements 
contribute to afjj (cf^), respectively. It then follows 
that the first term on the right hand side of Eq. 129|) 
corresponds to cr^ — c/Stt, whereas the second term 
(the vertex corrections) is the dissipative response with 

D e 1 N+T+Jl- N^T^J^V2-Va _ e 

(35) 

and only the sum of the reactive and dissipative response 
is zero. 

In summary, we calculated the spin-Hall conductiv- 
ity in a two dimensional electron gas within the self- 
consistent Born approximation, including the vertex cor- 



rections in the ladder approximation. We remark that, 
although a number of similar studies exist in the recent 
literature, the final conclusions are often contradictory. 
This may be due to the fact that the relevant integrals 
depend in a very subtle way on the type of the phys- 
ical limit considered. For this reason in this work we 
evaluated all the relevant integrals both analytically and 
numerically. This allowed us to confirm the conclusions 
of RefsAiiii^. In particular, we find that the spin-Hall 
conductivity is strongly suppressed below the universal 
value of e/87r. Furthermore we have demonstrated that 
the result is not only valid for pure s-wave scattering, 
but is robust upon the inclusion of a weak momentum 
dependence of the scattering probability. 
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